Indirect search for dark matter: prospects for GLAST 
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Possible indirect detection of neutralino, through its 7-ray annihilation product, by the forthcoming 
GLAST satellite from our galactic halo, M31, M87 and the dwarf galaxies Draco and Sagittarius 
is studied. 7-ray fluxes are evaluated for two representative energy thresholds, 0.1 GeV and 1.0 
GeV, at which the spatial resolution of GLAST varies considerably. Apart from dwarfs, which are 
described either by a modified Plummer profile or by a tidally-truncated King profile, fluxes are 
compared for halos with central cusps and cores. It is demonstrated that substructures, irrespective 
of their profiles, enhance the 7-ray emission only marginally. The expected 7-ray intensity above 1 
GeV at high galactic latitudes is consistent with the residual emission derived from EGRET data if 
the density profile has a central core and m x < 50 GeV, whereas for a central cusp only a substantial 
enhancement would explain the observations. From M31, the flux can be detected above 0.1 GeV 
and 1.0 GeV by GLAST only if m x < 300 GeV and if the density profile has a central cusp, case in 
which a significant boost in the 7-ray emission is produced by the central black hole. For Sagittarius, 
the flux above 0.1 GeV is detectable by GLAST provided the neutralino mass is below 50 GeV. From 
M87 and Draco the fluxes are always below the sensitivity limit of GLAST. 
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I. INTRODUCTION 

The nature of dark matter (DM) remains a major unre- 
solved problem in astroparticle physics. Main evidences 
for dark matter come from astrophysical observations. 
The rotation curve of a typical spiral galaxy flattens 
or slowly rises beyond the optical radius which includes 
most of the galaxy's luminosity [1,2]. The flattening of 
the rotation curves indicates the presence of DM at large 
radii (a few tens of kpc) : a result which is confirmed, up 
to scales of 30 kpc, by the study of close binary systems 
[3]. Among other evidences is the confinement of hot 
X-ray coronae around giant ellipticals which requires a 
massive dark halo [4,5]. At the centre of spirals or giant 
ellipticals, kinematical data indicate a large ratio of the 
baryon to dark matter density (B/DM) [6-8]. On the 
contrary, for dwarf galaxies it seems that the presence of 
a substantial amount of DM in the very central regions 
is required to explain the kinematical data [9]. On large 
scales, the average mass to light ratio (M/L) for galaxy 
clusters [10], the infall of the Local Group towards the 
Virgo cluster [11], the velocity dispersion of galaxies [12] 
and the angular power spectrum of the cosmic microwave 
background (CMB) radiation [13] have been long provid- 
ing evidences for dark matter. Various CMB experiments 
and distances to type la supernovae favor a flat universe 
of total density parameter f2^ ~ 1 and the matter density 
parameter il m w 0.30 [14,15]. The total baryonic con- 
tribution, fixed either by the primordial nucleosynthesis 
of light elements or by the secondary acoustic peaks ob- 
served by CMB experiments such as WMAP is f2b ~ 0.04 
[13], indicating that on large scales gravitational forces 
are mainly due to non-baryonic matter. In view of all 
these evidences, the inevitable question arises : what is 
dark matter made of? 

Most dark matter candidates which arise in the stan- 
dard model of particle physics are excluded by the 



present observational constraints. For instance, neutri- 
nos (v) which decouple relativistically from the primor- 
dial plasma and have a high relic abundance are the pri- 
mary hot dark matter (HDM) candidates. However, in a 
HDM Universe, small-scale structures, of typical galaxy 
sizes, are erased by relativistic streaming, hence neutri- 
nos are ruled out as DM candidates. Moreover, WMAP 
data set a robust upper limit $l v h? < 0.0076 to their den- 
sity [13]. Supersymmetric (SUSY) extensions of the stan- 
dard model lead to many more new DM candidates such 
as s-neutrinos, axions, gravitinos, photinos and cryptons 
[16]. Presently, the most plausible SUSY dark matter 
candidate is neutralino (x) which is the lightest super- 
symmetric particle. Neutralino is stable and hence is a 
candidate relic from Big Bang, if i?-parity quantum num- 
ber, introduced to avoid a too rapid decay of proton, is 
conserved as is the case in the Minimal Supersymmetric 
extension of the Standard Model (MSSM). Neutralino 
is an electrically neutral Majorana fcrmion whose mass, 
can range from a few GeV to a few hundreds of TcV. 

~ 30 GeV has been set by the 
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LEP accelerator [17], while an upper limit of m x ~ 340 
TeV is favored theoretically to preserve unitarity [18]. 

In the past few years, major experiments have gone 
underway for the direct or indirect detection of dark mat- 
ter particles [19]. Direct-detection experiments, such as 
DAMA, EDELWEISS, CRESST and IGEX are among 
many others [20] , which basically measure the energy (up 
to tens of keV) deposited at the detector by the elastic 
scattering of dark matter particles from the detector nu- 
clei. Direct detection experiments also use the annual 
modulation of the signal due to the orbital motion of the 
Earth around the Sun in their DM search. Clearly, the 
event rate, the search strategy and data analysis strongly 
depend on the spatial distribution of dark matter and 
its dynamics in the galactic halo, which are not unfor- 
tunately well-understood. For instance, it is not estab- 
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lished if dark matter halos are dynamically-relaxed struc- 
tures or not. Whether dark matter is homogeneously dis- 
tributed with isotropic velocity distribution or whether 
there are inhomogeneities such as local streams, e.g. like 
that manifested through the tidal arms of the Sagittar- 
ius dwarf [21,22] is not entirely clear. The leading trail of 
Sagittarius could be showering matter down upon the So- 
lar System, affecting the local halo density and perturb- 
ing the velocity distribution [23]. The presence of very 
high density structures such as caustics, if survived up to 
the present day, could also have dramatic effects on di- 
rect searches [24]. Moreover, dark halos are generally not 
at rest and have considerable angular momentum [25], 
whose vector direction is probably not the same as those 
of the present spin axes of (baryonic) disks [26] . If this 
was the case for the Galaxy, then significant variations 
in the modulation of the signal would be expected. All 
of these are just a few of many uncertainties about prop- 
erties of DM in the halos which overshadow the search 
experiments. 

The indirect detection experiments search for products 
of self-annihilation of neutralinos such as energetic lep- 
tons, hadrons and also particles which would emerge in 
the follow-up hadronization and fragmentation processes. 
In addition to 7-ray lines generated through the annihila- 
tion channels XX — > 77 an d XX — * the annihilation 
of two neutralinos also produces a 7-ray continuum as 
a consequence of the neutral pion decay. Besides these 
energetic 7s, neutrinos are also produced either in quark 
jets (bb interactions) or in the decay of t leptons and 
gauge bosons. Neutrinos produced in former process are 
less energetic than those produced in the latter. Neutrali- 
nos can be decelerated by scattering off nuclei and then 
accumulating at the centre of the Earth and or at the cen- 
tre of the Sun (or inside any other gravitational potential 
well), thus increasing the annihilation rate. However, so 
far different experiments designed to detect decay prod- 
ucts like high-energy neutrinos, have only managed to set 
upper limits on fluxes coming from the Earth centre or 
from the Sun [27]. 

The 7-ray emission on the other hand is expected, for 
instance, from the galactic centre and other massive dark 
matter halos. Since the annihilation rate depends on the 
square of the local density, the distribution of DM in 
the halo is of crucial importance. The degree of dumpi- 
ness, the density profile of a dark matter halo and in 
particular whether it has a central cusp or not, as well as 
the presence or absence of a central supermassive black 
hole (SMBH) could influence the annihilation rate. Other 
events such as the infall of small satellites, which are not 
totally disrupted by tidal forces can locally enhance the 
7-ray emission [28]. 

The prediction of 7-ray fluxes requires two separate in- 
puts: that coming from particle physics for issues such 
as the interaction cross section and the number of pho- 
tons per annihilation, and the input from astrophysics for 
problems such as the spatial distribution of dark matter 
in potential sources. 



Most studies on indirect DM detection explore various 
decay channels of neutralino annihilation in the huge pa- 
rameter space of MSSM, consisting of 91 real parameters 
and 74 phases (most of them can be absorbed by field 
redefinitions [29]). In the effective MSSM model, where 
the low-energy parameters are constrained by accelera- 
tor data, seven free parameters still remain (the higgsino 
mass parameter /i, the gaugino mass parameter M2, the 
ratio of the Higgs vacuum expectation values tan/3, the 
mass of the CP-odd Higgs boson tha, the scalar mass 
parameter m and the trilinear soft SUSY-breaking pa- 
rameters Af, and A t ). The computation of SUSY particle 
spectrum is therefore a difficult and sometimes uncertain 
procedure, which is often done using publicly available 
numerical codes like SUSPECT able to explore the full 
parameter space [33]. In this work, we take a different 
approach and assume that the number of photons pro- 
duced per annihilation, Q 7 , can be evaluated under the 
assumption that the XX annihilation process is similar to 
a QCD jet. In this picture, the evaluation of Q 7 is rather 
simple since it depends only on one scale parameter, the 
neutralino mass. 

The astrophysical aspect of the neutralino flux calcu- 
lation, which concerns the dark matter distribution, re- 
mains rather uncertain. Numerical simulations suggest 
that the density in the central regions of dark halos varies 
as p cx 1/r 13 , where (3 takes on values such as f3 = 1 [34] 
or [3 = 3/2 [36]. However, as we mentioned earlier, a cusp 
profile is not always supported by the observational data 
such as the rotation curves of bright galaxies [1] . More- 
over, in order to avoid the inevitable divergence of the 
7-ray emission rate calculated from these cusp profiles, 
various assumptions about the central density or cutoff 
radius are always made (see for example [28,37,38]). As 
a consequence, the value of the square of the density in- 
tegrated along the line of sight in the direction of the 
galactic centre performed by different authors varies by 
as much as four orders of magnitude (see Table 1). The 
density profile and the M/L ratio also depends on the 
type of the galaxy, e.g., as compared to normal spirals, 
dwarfs seem to have a much higher M/L ratio and do not 
seem to have a central cusp. Thus, as far as the selection 
of a suitable source is concerned, many factors such as the 
possible density profile of the source, the level of back- 
ground contamination, the distance to the source and the 
M/L ratio should be taken into account. Different studies 
have considered the galactic centre direction as a privi- 
leged source due to its high column density. However, the 
7-ray emission from our halo is highly contaminated by 
the local background, mostly produced by cosmic ray in- 
teractions with the interstellar environment. Background 
contamination is less significant for sources outside the 
galactic plane and extragalactic sources which subtend a 
small solid angle at the detector. Because of their prox- 
imity, the galactic centre and the nearby dwarfs, Draco 
and Sagittarius, offer unique sites for dark matter detec- 
tion. Due to their considerably large dark matter con- 
tent, M31 and in particular M87, are also potential 7-ray 



2 



sources. These are the sources that we shall consider in 
this article. We have calculated 7-ray fluxes from the 
Milky Way, M31 and M87 using two different density 
profiles: the first one is a very steep but non-diverging 
central profile based on very recent numerical simulations 
[39], while the second is a core-profile (Plummer profile) 
chosen because the energy distribution, the density and 
the potential are known in a closed form. The Plummer 
profile also leads to a finite mass and the density does not 
diverge at the centre. For Draco, we have used a density 
profile resulting from a dynamical model whereas a King 
profile was adopted for Sagittarius. In addition to these 
computations, we have simulated halos with clumps in 
order to verify how substructures affect the predicted 7- 
ray fluxes and have found that dumpiness enhances the 
predicted emission only marginally, in disagreement with 
some previous claims [37,66] and in agreement with oth- 
ers [38]. 

Unlike most previous studies which estimate fluxes 
mainly for atmospheric Chcrcnkov telescopes (ACTs) 
such as CELESTE and VERITAS, in this work, we con- 
sider the possibility of detection by the future Gamma- 
ray Large- Area Space Telescope (GLAST), which has 
various obvious advantages over ACTs, the more impor- 
tant of which are as follows: i) lower energy threshold, 
allowing to probe neutralino masses above 10 GeV; ii) 
the background is mainly due to the diffuse extragalac- 
tic emission and iii) the spatial resolution varies with the 
threshold energy, allowing to probe the halo density pro- 
file. 

Our results suggest that GLAST may detect M31 at 
energies above 0.1 GeV and 1.0 GeV if the neutralino 
mass is less than 300 GeV and only if the profile is cusped. 
Core profiles require neutralino masses lower then up- 
per limits established by accelerator data. Thus GLAST 
can put constraints not only on the neutralino mass but 
also on the density profile, combining the different spatial 
resolution at different energies. According to our compu- 
tations, in spite of having a very massive halo, emission 
from M87 cannot be detected either above 0.1 GeV or 1.0 
GeV. However, for M87, there are enough evidences for 
the presence of a central SMBH [40,41] which can form 
a central dark matter "spike" and which would conse- 
quently boost the 7-ray flux by a factor of about 200, 
producing a detectable signal. 

For the dwarf spheroids Draco and Sagittarius, we eval- 
uate the flux and demonstrate that 7-ray emission from 
Draco is below the sensibility limit of GLAST, whereas 
Sagittarius can be detected above 0.1 GeV energy thresh- 
old, if the neutralino mass is < 50 GeV. 

This article is organized as follows. In Section II, we 
discuss the 7-ray emission and the flux equation. The an- 
nihilation rate and the number of photons produced per 
annihilation event are discussed in Section III. In Section 
IV we discuss the astrophysical parameters required for 
the flux calculation. In Section V, we present our nu- 
merical results for emission from the galactic halo , M31, 
M87, Draco and Sagittarius. In Section VI, we demon- 



strate how dumpiness of halos would effect the 7-ray flux. 
In Section VII we study flux enhancement effects due to 
the central SMBH in M31 and M87 and finally in Section 
VIII, we summarize our results and conclude our work. 



II. THE 7-RAY FLUX 

Let us assume that dark matter halos are constituted 
by neutralinos: the lightest supersymmetric particle in 
many SUSY models (for a review of SUSY dark matter 
sec [42]). Neutralino is a linear superposition, 



(1) 



of four mass-eigenstates where B and W are gauginos 
(the superpartners of electroweak gauge bosons) and 
arc higgsinos. If the sum a 2 - + a 2 ~ 7 > 0.9, the neutralino 
is gaugino type, while it is higgsino type if that sum is 
less than 0.1. 

We further assume that in the early universe, neu- 
tralinos were in thermal equilibrium with the primordial 
plasma and that their present density was fixed at the 
freeze out, i.e., at the temperature T* when they decou- 
pled from the other particle species. At freeze out, the an- 
nihilation rate is comparable to the expansion rate of the 
Universe and moreover neutralinos are non-relativistic 
(kT* < m x c 2 ). After decoupling, neutralinos can an- 
nihilate via the following channels 

XX -» IT, qq, W + W-,Z Z°,H°H , Z°H°, W ± H^ (2) 

into leptons and anti-leptons (11), quarks and anti-quarks 
(qq), charged (W ± and and neutral (H° and Z°) 

bosons respectively. 

The decay of neutral pions formed in the hadroniza- 
tion process is the dominant source of continuum 7-rays. 
Besides the continuum emission, two annihilation chan- 
nels may produce 7-ray lines. The first, XX ~> 77; where 
the photon energy is <~ m x and the second, XX ~^ 
where the photon energy satisfies e 7 = m x — m|/(4m x ). 
The latter process is only important for neutralino masses 
higher than ~ 45 GeV. 

Since neutralinos are Majorana particles, their density 
is equal to that of anti-neutralinos and hence the annihi- 
lation rate per unit volume is 



Px 



(3) 



where < <J xx v > is the thermally- averaged annihilation 
rate, m x is the neutralino mass and p x is the neutralino 
matter density. If the source is spatially extended, the 
radiation intensity Ij(r p ), measured at a given projected 
distance r p from the centre is 



(4) 
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where 



1 = 



(5) 



is the reduced intensity, Q 7 is the number of photons 
above a certain energy threshold produced per annihi- 
lation and the integral is carried along the line of sight. 
The flux / 7 within a solid angle AQ subtended by the de- 
tector (dfl — 2TTr p dr p /D 2 for a small angular distance) 



is 



f-y = J I-ydfl = Jj2 J Iy(r P )P(r p )r p dr p 



(6) 



where D is the distance to the source and P{r p ) is the 
point-spread function (PSF) of the detector. In the case 
of GLAST, we adopt a Gaussian PSF (P(r p ) = e^l' 2 * 2 ) 
whose width depends on the energy. The GLAST project 
specification gives the width at which 68% of the sig- 
nal is included, which yields a(0.1GeV) w 1.05° and 
a(lGeV) w 0.17°. 



III. NEUTRALINO ANNIHILATION RATE : THE 
VALUES OF < <j xx V > AND Q 7 

In order to compute the 7-ray flux from (6), different 
physical and astrophysical parameters should be calcu- 
lated. First, the average annihilation rate, which will be 
estimated under the assumptions already mentioned, i.e., 
that neutralinos are initially in thermal equilibrium and 
decouple at T* , when the annihilation rate equals the ex- 
pansion rate. The latter is determined by the Friedman 
equation relating the Hubble parameter to the total en- 
ergy density. As we shall see, neutralinos decouple at 
temperatures in the range 0.4 — 70 GcV, depending on 
their assumed mass. Once neutralinos decouple, their co- 
moving number remains almost constant, but their con- 
centration with respect to photons changes, since other 
particles which decouple later, produce a "reheating" of 
the photons due to entropy conservation. In this case, the 
ratio between the temperature after and before a given 
decoupling depends on the variation of the effective num- 
ber of degrees of freedom as (g e f f,bef / 9eff,aft) 1 ^ 3 , where 



9eff(T) = ±J2(g B + -g F ) 



(7) 



and the sum is over all bosons, B, and fermions, F, 
present in the primordial plasma. For temperatures 
above 0.4 GeV and below 80 GeV one would expect 
that quarks and gluons arc present and that the gauge 
bosons (W^ and Z°) have already disappeared. In this 
case, the main particles contributing to g e ff are pho- 
tons, gluons, the three leptons and their anti-particles 
(the tauon decouples at T* <~ 1.8 GeV) as well as their 
associated neutrinos and anti-neutrinos, the six quarks 
and their anti-quarks. 



Under these conditions one obtains (see for instance 
reference [43]) the balance equations 



.. u2 6.88 x 10 7 T . , /2 _ v 

fi x>4 = 7 7^m x (GeV)y 3 ^e » , 



9eff(T*) 



(8) 



8.82 x io- 3 y ge// (r.) _ 1/2 3 _! 

dvy!) = T^r-rrs V e cm s , (9) 

V xx ' g x m x (GeV) w 

where y = m x c 2 /kT* and h 65 is the Hubble parameter 
in units of 65 kms _1 Mpc _1 . From these equations, once 
the density parameter is fixed, the decoupling tempera- 
ture and the annihilation rate per particle can be com- 
puted for a given ncutralino mass. Here we adopt the 
cosmologically accepted value of £l x — 0.26 and assume 
that the neutralino mass lies in the range 10 GcV< m x < 
2000 GeV. In this case, the decoupling temperature pa- 
rameter y varies within the interval 22.5 < y < 28.2, cor- 
responding to the temperature range mentioned above. 
The mean annihilation rate per particle varies very little, 
(7.7 — 9.5) x 10 _27 cm 3 s _1 and is comparable to values 
adopted in other studies [28,44]. 

A major difficulty in the calculation of the number 
of 7s per annihilation stems from the unclear details of 
the hadronization and gamma-ray emission in the large- 
dimensional parameter space of SUSY models. Here, we 
use an approximation to the numerically computed QCD 
fragmentation functions taken from [45,46]. We assume 
that hadronization and subsequent 7-ray emission is sim- 
ilar to the decay of QCD jets and that the annihilation of 
two neutralinos results in two jets, each of energy ~ m x . 
Using the fragmentation function of jets into photons, 
neutrinos and baryons, which is known to a good ap- 
proximation up to a few TeV [45,47], the number of pho- 
tons per annihilation and per jet produced in the energy 
interval x, x + dx is 



dQ 



dx 8 V 3 V ^ + 3 



(10) 



where the dimcnsionless energy is defined as x — 
e 7 /(m x ), with e 7 being the photon energy. Therefore, 
the total number of photons per annihilation event above 
the threshold energy e 7it h is simply given by 



1 d^y 

dx 



dx . 



(11) 



where the factor two takes into account the formation of 
two jets rather than one, due to momentum conservation. 

In Fig. 1, we show the photon production rate < 
<J xx v > Q~ ( in cm 3 s _1 for two threshold energies (0.1 
and 1.0 GeV) as a function of the neutralino mass. The 
rates, in spite of our simplifying assumptions, are in good 
agreement with estimates derived from SUSY numerical 
codes [28,44]. This is not surprising since the parameters 
required in SUSY codes are chosen to match the expected 
relic density. We have taken a reversed approach, i.e., we 
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fix the relic density first and then compute the annihi- 
lation rate using the balance equations (8) and (9). As 
a consequence, both procedures should give comparable 
results, but our method has the advantage of using only 
one free parameter, the neutralino mass. 



given by the expression (5) estimated by several authors. 
In fact, these reduced intensities correspond to values 
averaged over a given solid angle representative of the 
spatial resolution of the detector. 
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FIG. 1. The variation with the neutralino mass of the pho- 
ton production rate per annihilation above the energy thresh- 
olds, e 7 , th , 0.1 GeV and 1.0 GeV. 



IV. ASTROPHYSICAL PARAMETERS: THE 
VALUES OF T 



In the previous section, we have evaluated the annihi- 
lation rate < o~ xx v > an d the total number of photons 
per annihilation event above a given energy threshold, 
Qj. In this section, we evaluate the reduced intensity 1, 
given by (5) whose value is needed for the evaluation of 
the 7-ray flux (6). The integral in (5) is often evaluated 
assuming a density profile of the general form 



p(r) = 



(r/r*)P(l + r/r*p 



(12) 



where p* and r* are the characteristic density and ra- 
dius. If (3 = 1 and 7 = 2, one obtains the so-called 
Navarro- Frenk and White (NFW) profile [34], while if 
(3 = 7 = 1.5 one obtains the so-called Moore profile [36]. 
Clearly, when integrated along the line of sight passing 
through the centre, as required by expression (5), both 
of these two profiles diverge. Various tricks are used to 
overcome this difficulty. One can impose a finite central 
density such as p x — m x H n / < o xx v >, since denser 
regions will be depleted by high neutralino annihilation 
rates [37]. High resolution simulations are also used to 
study the behaviour of the density profile in the very 
central regions of halos. However, even the highest res- 
olution simulations are presently unable to resolve the 
very inner structure of halos and different extrapolations 
are needed to obtain the central density. In order to il- 
lustrate this point, Table 1 compares different values of 
density-squared along the direction of the galactic centre 



Table 1. Reduced intensity in the direction of the galactic 
centre 



Profile 


Jp 2 ds (GeV 2 cm" 5 ) 


Reference 


Moore 


3.3 x 10 26 


[37] 


NFW 


2.8 x 10 25 


[57] 


core 


3.0 x 10 22 


[57] 


cusp 


2.4 x 10 22 


[48] 


NFW 


5.2 x 10 25 


[38] 


SWTS 


1.8 x 10 24 


[38] 



Simple inspection of Table 1 shows that the integral (5) 
may vary by four orders of magnitude according to the 
assumed profile. Consequently, severe uncertainties over- 
shadow the value of the predicted 7-ray flux. 

Rather recently, the inner structure of dark halos in 
ACDM cosmology was revisited [39] and the new "a- 
profile" 



p{r) = p*exp 



1 



(13) 



fitting the the numerical results was proposed, which has 
several advantages over (12). Firstly, the total mass and 
the central density are finite which solve the aforemen- 
tioned divergence problem in the evaluation of the cen- 
tral intensity. Secondly, the logarithmic slope decreases 
inwards more gradually than NFW or Moore profiles. Ac- 
cording to simulations, the parameter a is restricted to 
be in the range 0.1 — 0.2 and here we adopt for M31 and 
the Milky Way the suggested value of a = 0.17 [39]. The 
two remaining parameters, and r*, can be estimated if 
the total halo mass and one other physical quantity such 
as the maximum circular velocity are known. However, 
the rotation curve within the effective radius is not dom- 
inated by dark matter [25] and so the maximum rotation 
velocity is not a robust quantity. The halo of M31, as 
we shall see, is probably similar to our own halo which 
implies that the dark matter energy density at a distance 
of about 8 — 9 kpc from the centre should be comparable 
to that observed in the solar vicinity, i.e. 0.4 GeVcm~ 3 . 
We use this value together with the total mass of M31 to 
fix the parameters p* and r* . 

Present observational data constraining the mass of the 
halo of M31 are rather poor. For instance, the rotation 
curve seems to be well-explained with just the baryonic 
matter in the bulge and the disk, whose mass distribu- 
tions are traced by the blue light [49]. This does not 
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exclude the existence of a dark matter halo whose con- 
tribution to the gravitational force, and hence its effect 
on the rotation curve, would be more relevant only at 
distances larger than few tens of kpc. Direct estimates 
of the halo mass of M31 are based on the peculiar mo- 
tion of the Local Group (LG) and on the kinematics of 
its satellites. A new solution for the solar motion which 
includes the 32 probable members of the Local Group 
and assumes virial equilibrium, leads to a total mass for 
the LG of (2.3 ±0.6) x 1O 12 M [50]. Using the projected 
mass method ( [51,52]), the contribution of the M31 sub- 
group (including only 7 satellites) to the total mass of 
the LG amounts to (13.3 ± 1.8) x 10 11 M [50]. A similar 
result (12.3tg 8 x lO n M ) was obtained [53], including 
10 satellites, 17 distant globular clusters and 9 halo plan- 
etary nebulae, while a firm upper bound of 24 x 10 11 M 
was derived from a similar analysis including 15 satellites 
and new radial velocity data [54]. Using the kinematical 
data on these 15 satellites and the projected mass estima- 
tor [52], we have obtained a mass of (1.5 ±0.5) x 1O 12 M 
for M31. Other indirect estimates based on disk forma- 
tion and cusp dark matter halo expected for ACDM cos- 
mology predict a mass of 1.6 x 10 12 M [55]. We adopt 
the value 1.5 x 1O 12 M , which is compatible with all of 
these different estimates. 

Subsequently, using the mass of M31 and also the local 
density of 0.4 GcVcm -3 which we argued for earlier in 
this section, we obtain the two parameters defining the 
"a-profile" (13): p* = 4.5 x lO^gcm" 3 and r* = 11.6 
kpc. We emphasize that these parameters give, in terms 
of the total potential energy | W \ , a gravitational radius 
of r g — GM 2 / | W |s=s 250 kpc, which is a typical value 
that we have found in our numerical simulations for halos 
of similar masses [25]. 

Since the high B/DM ratios observed in the central 
regions of galaxies is contradictory with a cusp in the 
DM density profile, we have also considered a profile with 
a central core to represent the dark matter distribution 
inside halos. We assumed a Plummer density profile, 

P[r) = [l + (l/3Kr/ro) 2 p/2' (W) 

and since the central density and total mass are finite, the 
distribution function depends only on the total energy 
( f(E) — K | E I 7 / 2 ) and the gravitational potential 
can also be obtained in a closed form. Using the two 
constraints mentioned earlier in this section, we obtain 
the two parameters defining the profile (14) in the case 
of M31: po = 2.8 x lO^gcm" 3 and r = 12.2 kpc. 

For M87, the properties of the gaseous X-ray corona 
were determined in reference [56] . From the density and 
temperature profiles of the hot gas, the total mass in- 
side a given radius r can be computed under the as- 
sumption of hydrostatic equilibrium. After correction 
for the contribution of baryons (stars and gas) to the 
gravitational forces, the dark matter distribution can be 
evaluated [56]. Remarkably, such a distribution can be 
well-represented by an "a-profile" with the parameters 



a = 0.20, p* = 1.07 x 10- 25 gcm- 3 and r* = 79.9 kpc. 
Integration of this profile gives a mass of 9.7 x 10 12 M 
inside a radius of 100 kpc, in agreement with other esti- 
mates [56]. However, we cannot exclude that at distances 
less than a few kpc, such a profile would be inconsistent 
with observations. 



Table 2. Values of the parameters used in the density profiles. 



M31 

cusp : a = 0.17 p* = 4.5 x 10~ 25 g/cm 3 
r» = 11.6 kpc 

core : p = 2.8 x 10~ 24 g/cm 3 

r = 12.2 kpc 

M87 

cusp : a = 0.20 p* = 1.07 x 10~ 25 g/cm 3 
r* = 79.9 kpc 

Draco 

Modified Plummer profile : 

= -0.34 
7 = -1.03 

po = 3.8 x 10~ 23 g/cm 3 

Sagittarius 

King profile : 

r t = 2 kpc 

r c = 0.55 kpc 

po = 2.2 x 10~ 24 g/cm 3 

(15) 

These parameters, summarized in Table 2 (for dSph's 
see details in Section V), will be used in our predictions 
of the expected 7-ray flux from neutralino annihilation in 
M31 and M87. In Fig. 2 we show the expected intensity 
profile for M31 and M87 as a function of the angular dis- 
tance to the centre. These profiles were derived numer- 
ically from (4), for an energy threshold of 0.1 GeV and 
m x = 50 GeV. The comparison uses only density profiles 
with central cusps, defined by the parameters discussed 
above. Note that although M31 has a less massive halo, 
its central brightness is expected to be higher than that 
of M87. 
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3 gives the expected intensity function of the neu- 
tralino mass for the a-profile (with a central cusp) and 
the Plummcr profile (with a central core). 

Table 3. Intensity (in photons cm~ 2 s" 1 sr^ 1 ) above 1 GeV in 
the direction b = 90° 



m x (GeV) 


a— profile 


Plummcr profile 


30 


1.1 x 10~ 8 


1.9 x 10~ 7 


50 


6.4 x 10~ 9 


1.1 x 10~ 7 


100 


2.8 x 10~ 9 


4.9 x 10~ 8 


200 


1.1 x 10~ 9 


2.1 x 10~ 8 


300 


6.8 x 10~ 10 


1.2 x 10~ 8 


500 


3.4 x 10~ 10 


6.1 x 10~ 9 



FIG. 2. Intensity profiles for M31 and M87 for photon en- 
ergies above 0.1 GeV. Intensities were computed using the 
cusp density profile (13) and m x = 50 GeV 

V. NUMERICAL RESULTS 
A. Flux from the galactic halo 

It is instructive to compare our results with previ- 
ous estimates of the reduced intensity (5) in the di- 
rection of the galactic centre. In our work, the inte- 
gral of the density squared along the line of sight is 
1.6 x 10 26 GeV 2 cm- 5 if the a-profile (13) is adopted. 
This compares quite well with the value of [37], based 
on extrapolation from high resolution simulations, but 
is one order of magnitude higher than those given in a 
few of the previous works [38,57], where NFW profiles 
were adopted. Using the core profile, this integral (5) is 
equal to 1.2 x 10 23 GeV 2 cm~ 5 , considerably smaller than 
estimates for profiles with central cusps as would be ex- 
pected, but still one order of magnitude higher than some 
others which also consider a core profile (see for example 
[57] whose estimate of X is given in Table 1). 

In the energy range 0.1 — 1.0 GeV, cosmic ray in- 
teractions give a substantial contribution to the 7-ray 
emission originated from the galactic disk. The electron 
component produces high energy photons either by in- 
verse Compton scattering or bremsstrahlung, while pro- 
tons produce 7-photons via the decay of neutral pions 
generated in collisions with interstellar matter. In a pre- 
vious work which uses EGRET data, at high galactic 
latitudes a significant residual intensity of 10 -7 — 10~ 6 
photons cm _2 s _1 sr _1 above 1 GeV is claimed even after 
correction for the expected background of cosmic rays 
and the diffuse extragalactic emission [58]. These resid- 
ual intensities although marginal have rather interesting 
implications for the halo density profile. We computed 
the expected intensity, averaged over cells of 0.5° x 0.5° 
as considered in [58], in the direction b = 90°. Table 



The results summarized in Table 3 show that at high lat- 
itudes, contrary to what is obtained in the direction of 
the galactic centre, the Plummcr profile produces an in- 
tensity higher than that derived from the a-profile. This 
is easily understood since the latter profile gives a larger 
mass concentration near the centre while the former has a 
shallower mass distribution. Furthermore, Table 3 shows 
that intensities derived from the a-profile are always be- 
low the EGRET residual intensity. In this case, an im- 
portant intensity enhancement by clumps is necessary to 
explain the EGRET residuals [66,37]. Intensities derived 
from the Plummer profile are comparable to EGRET val- 
ues if m x < 50 GeV, a limit not inconsistent with the 
lower bound derived from LEP. As we shall see in Sec- 
tion VI, the enhancement of the 7-ray emission produced 
by subhalos is rather small, with boost not exceeding a 
factor of 2. Presently, a firm conclusion cannot be made 
since the EGRET residuals are in the sensibility limit of 
the instrument. This situation is expected to improve 
greatly with the forthcoming GLAST. 

B. Fluxes from M31 and M87 

The detection of 7-rays with energies around 50 GeV 
from M31 was considered by the CELESTE collabora- 
tion [30] . Their study shows that the annihilation signal 
could be detectable by CELESTE for m x > 200 GeV, if 
the M31 halo is very clumpy and/or if there is a strong 
accretion onto a central supermassive black hole [30] . 

M87, the giant elliptical galaxy at the centre of the 
Virgo cluster, is also a potential source of 7-rays pro- 
duced by XX annihilation [31]. If it is simply assumed 
that the neutralino mass is 1 TeV and that the photon 
production rate above the threshold energy of 50 GeV 
is < <J xx v > Q 7 = 1.5 x 10 _24 cm 3 s _1 , the predicted 7- 
ray flux from M87 would be below the sensibility of the 
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present ACTs [31]. However, detection by the new gen- 
eration of ACTs like VERITAS would be possible if there 
is an enhancement by a factor 40 due to the dumpiness 
of the halo [31]. In this work, we consider the possibil- 
ity of detection by GLAST, which has various obvious 
advantages over the atmospheric Cherenkov telescopes, 
which we have mentioned in Section I. 

The photon production rate above a given energy is 



i? 7 =< cr xx v > Q~f = 



47T/. 



7, min 



/ P(r p )dCl J p%ds 



(16) 



where / 7 , m in is the minimum flux above a certain photon 
energy detectable by GLAST and we have used equations 
4 and 6. Values of / 7 , m in varies from one 7-ray experi- 
ment to another and here, we use the sensibility curves 
summarized in [59] . The production rate i? 7 given by the 
above expression represents the minimum required value 
for the production of a detectable signal. In Fig. 3 we 
plot the minimum photon production rate and the pho- 
ton production rate derived in Section III as a function of 
the neutralino mass, for an energy threshold of 0.1 GeV. 
The required minimum production rates are shown for 
both the a-profile and the Plummcr profile. Numerical 
calculations were performed adopting a distance of 770 
kpc for M31. 
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FIG. 3. Minimum photon production rate required for de- 
tection and theoretical estimates as a function of the neu- 
tralino mass and for an energy threshold of 0.1 GeV. Min- 
imum production rates are shown for the halo of M31 with 
cusp (a-profile ) and core (Plummer) profiles. 

Fig. 3 demonstrates that detectable fluxes from M31 
can be produced by both density profiles if m x < 25 GeV. 
In fact, the relatively large acceptance angle of GLAST 
at 100 MeV, includes a substantial fraction of the halo, 
reducing the expected differences due to different den- 
sity profiles. The situation is rather different when the 
threshold energy is 1.0 GeV (Fig. 4). In this case, the de- 
tector acceptance angle is considerably reduced and only 
a central cusp would lead to enough projected mass to 
produce a detectable flux if m x < 20 GeV. These mass 



limits are slightly smaller than the lower bound derived 
from LEP [17], but masses as low as 6 GeV have been 
derived from WMAP data in the frame of an effective 
MSSM model without gaugino-mass unification at the 
GUT scale [32]. In Section VII we discuss how these 
limits are changed by the presence of a central SMBH. 

In spite of its extremely massive halo, the predicted 
fluxes from M87 are about one order of magnitude below 
the sensibility limit of GLAST in the range 0.1-1.0 GeV, 
due to its large distance (17 Mpc) from us and in spite 
of the fact that a substantial fraction of a far halo would 
lie within the acceptance angle of the detector. As in the 
case of M31, the central SMBH strengthens considerably 
the emission, producing a detectable signal, as we shall 
shortly demonstrate in Section VII. 
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FIG. 4. Minimum photon production rate required for de- 
tection and theoretical estimates as a function of the neu- 
tralino mass and for an energy threshold of 1.0 GeV. Min- 
imum production rates are shown for the halo of M31 with 
cusp (a-profile) and core (Plummer) profiles. 



C. Fluxes from Dwarf Spheroidals 

Accumulating radial velocity data on stars of dwarf 
spheroidal (dSph) in the past years, have been providing 
convincing support for a large M/L ratios and hence a 
substantial amount of dark matter in these galaxies. The 
possibility of detecting 7-ray emission from dSph's due 
to neutralino annihilation has been considered in some 
recent works [31,57], most of which emphasize on the en- 
ergy range covered by ACTs. In the photon energy range 
of around 0.1 GeV, the acceptance angle of GLAST cov- 
ers a considerable fraction of the galaxy, hence increasing 
the amplitude of the received signal. 

Once again, the predicted flux depends on the dark 
matter density profile of dSph's. The light distribution 
in these systems is well-described by isotropic King mod- 
els [62]. However, in the case of Draco the velocity dis- 
persion rises slightly with the increasing radius, ruling 
out the possibility that the DM distribution be similar 
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to that of light [60]. The observed velocity dispersion 
profile can be reproduced by models in which the dark 
matter distribution is significantly more extended than 
light, having also a tangential anisotropy in the velocity 
dispersion [60]. These models are based on a modified 
Plummer potential and density family, characterized es- 
sentially by three parameters: ro,/3, and 7. The first de- 
fines the core or the scale radius, the second determines 
the mass distribution (f3 = 1 implies that DM follows 
light) and the third measures the velocity anisotropy. 
The potential is given by [61] 

0M = <Mi + i^ /2 (17) 

where (j>o is related to the central matter density by the 
equation po = 3 | (3 | 4>$ / [AnGr"^) . The central density 
can also be explicited in terms of the projected central 
velocity dispersion a as [61] 

2 

po = ^^m^) as) 

where F(j3, 7) = (5 + j3 - i)T{h/2 + P/2)/T{2 + (3/2). 

According to [60], the best parameters describing the 
Draco system are: r = 9'. 6 (0.2 kpc at a distance of 72 
kpc), /3 — -0.34 and 7 = -1.03 (implying tangential ve- 
locity anisotropy). Since the observed central projected 
velocity dispersion is 8.5 km.s -1 , from (18) one obtains 
for the central density p = 3.8 x 10~ 23 g.cm~ 3 . 

We have used this modified Plummer model and the 
parameters above to estimate the 7-ray flux from Draco. 
If we adopt m x — 30 GeV, corresponding to the lower 
limit derived from accelerator data, the expected 7- 
ray flux from Draco above 0.1 GeV is about 1.3xl0~ 10 
cm _2 .s _1 , which is about one order of magnitude less 
than the sensibility of GLAST at this energy threshold. 
Higher neutralino masses will give still lower fluxes. A 
similar result is obtained for 1.0 GeV energy threshold. 

The situation is more complicated in the case of Sagit- 
tarius since this dSph is partially disrupted by the tidal 
forces. Up to now, the different proposed models are not 
able to reproduce adequately the age and the observed 
structure of this dSph. Most of the models suggest that 
the system is disrupted after 1-2 orbits while observations 
indicate ten or more. 

Numerical simulations aiming to describe the merger 
history of this dwarf assumed either a Plummer [22] or a 
King profile [64] for the progenitor. These studies suggest 
that the mass of the core is in the range (5 — 20) x 10 s M Q 
[63,22,64]. 

Here we model the mass distribution of Sagittarius by 
an isotropic King profile, since no significant velocity dis- 
persion gradient is detected across the main body of the 
galaxy. The tidal radius is ill-defined in this system, but 
it should be smaller than 4 kpc in order to account for 
the 2 — 3 Gyr old M giants which were formed in the 
core and then escaped through the tidal boundary. If we 



take respectively for the tidal and core radii r t — 2 kpc 
and r c = 0.55 kpc, the total mass is about 1.8 x 10 9 M Q , 
consistent with simulations by [64]. Since the central 
projected velocity dispersion is 11.4 km.s -1 [63], the re- 
sulting central density is p n = 0.03 M Q pc~ 3 . In Fig. 5 
we show the minimum photon production rate required 
for detection and the theoretical rate as a function of the 
neutralino mass. Emission from Sagittarius can be de- 
tected by GLAST above 0.1 GeV, if the neutralino mass 
is < 50 GeV, but not above 1 GeV, due to the small 
acceptance angle at these energies. 
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FIG. 5. Minimum photon production rates and theoretical 
rates for energies above 0.1 GeV and 1.0 as a function of the 
neutralino mass. 



VI. EFFECTS OF CLUMPS ON THE FLUX 

Numerical N-body simulations show that dark matter 
halos contain a large number of self-bound substructures, 
which correspond to about 5 — 15% of their total mass 
[35,65]. These substructures are the consequence of the 
capture of small satellites, which have not yet been dis- 
rupted by tidal forces. Since the 7-ray emissivity depends 
on the square of the density, these substructures enhance 
the 7-ray flux relative to that expected for a smooth halo. 
However, whether this enhancement is significant or not 
is still a matter of debate. In some simulations, an im- 
portant enhancement of 7-ray flux due to dumpiness was 
reported for emission from our halo [37], while other nu- 
merical simulations seem to have indicated only a very 
small boost of the signal [38]. The enhancement due to 
clumps not only depends on the density profile of the 
clumps but also on their number density. The predicted 
number of massive substructures exceed the number of 
observed satellites of the Milky Way or M31, by at least 
one order of magnitude, which has been causing severe 
problems for ACDM model. For instance, a number of 
subhalos with masses above 10 8 M as high as 500 has 
been assumed in recent works [28]. 
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Here, to study the effects of dumpiness and the ef- 
fects of the density profile of clumps on the flux, we have 
made many different experiments with many test halos 
that we have produced numerically. In this way, we have 
been able to study the effects of dumpiness for halos with 
very high resolution and with various different profiles. 
Our test halos have the characteristic parameters of M31, 
with typical mass of Mh — 1.5 x 10 12 M Q , constituted 
by 5 x 10 7 particles, corresponding to a mass resolution 
of 3 x 10 4 M©. In the first experiment, the particles are 
distributed in a spherical volume of radius 250 kpc, com- 
parable to the gravitational radius expected for halos of 
such a mass (see [25]), according to a-profile. This will 
be our reference halo. In the next trials, we consider ha- 
los with substructures. The number of clumps N c j in the 
mass range m - m+dm is assumed to obey [35,38] 



dN a 



m 



1.78 



dm . 



(19) 



The normalization constant A is calculated by requir- 
ing the total mass in the clumps to be 10% of the halo 
mass and by assuming subhalos masses in the range 10 5 5 



10' 



9.25 



M Q . The minimum number of particles in the 



clumps is about 14, while the more massive ones have 
45, 000 particles. These figures correspond to mass bins 
built as described below. In spite of the low number of 
particles in the small clumps, our results are not seri- 
ously affected, since most of the flux enhancement is due 
to the massive structures. For numerical convenience, 
the number of clumps was calculated from (19) within 
logarithmic bins of width equal to 0.25. 

Clumps were distributed according to the (normalized) 
probability distribution p(r)d 3 r = (p(r)/Mh)d 3 r, where 
p(r) is assumed to have an a-profile, which gives the 
probability to find a clump at a distance r within the 
volume element d 3 r. The stripping process caused by 
tidal forces seems to reduce the density of a clump at all 
radii and, in particular, in the central regions, producing 
a density profile with a central core [67] . This result was 
further confirmed by simulations which found that the in- 
ner structure of subhalos are better described by density 
profiles shallower than NFW [38]. However, other simu- 
lations seem to indicate that the central regions of clumps 
are well-represented by power law density profiles, which 
remain unmodified even after important tidal stripping 
[68]. To test the importance of these differences for the 
flux, we allow for both of these possibilities: we simu- 
late both test halos with clumps having a central core 
p{r) oc l/(ro + r) 2 or NFW profile and study how our 
results would change from one to the other. 
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FIG. 6. Two-dimensional thin slices of the halos contain- 
ing 5xl0 7 particles, which we have numerically modeled to 
study effects of dumpiness on the 7-ray emission. The top 
panel shows one of our reference halos, which has the charac- 
teristics of M31, with an a-profile and no substructures. The 
bottom slice shows a halo with 10% of its mass distributed 
in substructures which, in the particular 2-dimensional plot 
show here, have NFW profiles. The number of substructures 
is taken from the power-law (19) and their centres are dis- 
tributed according to an a-profile. The clumps are highlighted 
by reducing the number of points of the background halo par- 
ticles. 

We recall that the 7-ray intensity (4) depends on the 
the integral along the line of sight of the square of the 
density. We define the enhancement factor to be the ratio 
between the reduced intensity calculated for halos with 
substructures to that calculated for our reference halo 
(with no substructure).. The numerical evaluation of in- 
tensities was performed by dividing the halo volume in 
a large number of small cubes. The cube volume varies 
logarithmically, being smaller in the central regions in or- 
der to increase the accuracy in the particle density. For 
angular distances 9 > 1°, the number of cubes is about 
10 6 whereas near the centre (9 < 1°) this number in- 
creases up to 15 x 10 6 . We have checked our algorithm 
for computing the reduced intensity X by comparing from 
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direct integration of (13). Errors are of the order of 2-3 
% if the integration is carried along line of sights far from 
the centre, but can reach values of up to 10% if the inte- 
gration is performed along directions close to the centre, 
due to resolution problems. Intensities were calculated 
in steps of A0 = 0.1° and Aip = 1.0°, with and ip being 
respectively the angular distance to the centre and the 
azimuthal angle around the centre of the galaxy. 
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FIG. 7. Plots show the variation of the enhancement factor 
as a function of the angular distance (in degrees) from the 
centre, for halos with characteristics of M31. Clumps have 
either an isothermal (top plot) or a NFW profile (bottom 
plot). In each case, three different simulations are shown to 
illustrate statistical fluctuations. 

In figure 7 we show the maximum enhancement factor 
as a function of the angular distance to the centre of M31 
for some selected trials. In the upper panel clumps have 
profiles with central cores whereas in the lower panel a 
profile with central cusp was adopted. We notice that, 
up to distances of about 4° — 5°, higher than the reso- 
lution of GLAST, independent of the clump profile, the 
enhancement is very small (about 10%), in agreement 
with the previous results (e.g. see [38]). Larger enhance- 
ments may be obtained at the outskirts of the halo and, in 



this case, clumps with central cusps produce a strength- 
ening almost one order of magnitude higher than clumps 
with cores. However, two points must be emphasized: a) 
these enhancements are local and when convolved with 
the PSF of the detector, they will be attenuated; b) as 
shown in Fig. 2, at 1° from the centre the intensity has 
already decreased by almost four orders of magnitude. 
Thus, at the halo boundaries, even with an enhancement 
factor of 40 — 50, the expected flux received by a de- 
tector pointing away from the centre will be below the 
dctcctability limit. 



VII. EFFECTS OF CENTRAL SUPERMASSIVE 
BLACK HOLES 

Fast-growing observational evidences indicate that 
most elliptical galaxies and bulges of spirals harbor 
SMBHs at their centres. In particular, recent spectro- 
scopic data suggest the existence of SMBH in M31, M87 
and in the Milky Way [70] . These kinematical data, if in- 
terpreted in terms of Keplerian motions, indicate a mass 
of 4.5 x 10 7 Mq for the SMBH in the centre of M31 [70] 
and of about 2.8 x 10 9 M Q for that in the centre of M87 
[40,41]. Kinematical data on stars very near Sagittarius 
A, the radio source supposed to be in the galactic centre, 
represent the best evidence in favor of the existence of a 
SMBH in the Milky Way [71], whose mass is 2.6 x 10 6 
M . 

If a black hole grows inside a dark matter halo, sup- 
posed to be constituted by a dissipationless fluid, the 
density around it also grows since the gravitational at- 
traction of the black hole causes a shrinking of the orbits 
inside the sphere of influence, whose radius is given by 
fbh = GMbh/c 2 , where a is the initial one-dimensional 
velocity dispersion of dark matter. If the growth is adia- 
batic [72], that is if tp << Mbh/Mbh << £r, where tp is 
the orbital timescale, Mbh the black hole accretion rate 
and £r the relaxation timescale, then a central spike will 
be formed inside the sphere of influence, having a central 
cusp profile of the form 



A 



(20) 



The value of the exponent 7 depends on the initial state 
of the matter distribution. On the one hand, if the ini- 
tial density profile has a central core, then the resulting 
spike is a power law with 7 = 3/2 [74,72]. On the other 
hand, an initial power law profile with an exponent 70 
will produce a steeper density profile with an exponent 
7 = 2 + 1/(4-70) [75]. 

The adiabatic scenario is questionable: firstly because 
the predicted accretion rates are considerably higher than 
the current estimates [73] and secondly, because mas- 
sive galaxies like M31, M87 and the Milky Way have 
most probably experienced several merging episodes [25] , 
which may have perturbed or even destroyed the central 
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spike. Nevertheless, we will assume in our calculations 
that there is a putative spike formed around the SMBH 
under adiabatic conditions. 

For M31, previous simulations [30] gives a spike profile 
p sp oc r~ 1,5 which is formed in less than 10 6 yrs. Ac- 
cording to [30] this profile gives a small enhancement of 
about 45%, with respect to a NFW profile, of the 7-ray 
emission from M31. Here we adopt the following proce- 
dure. We assume a smooth transition at rbh from the 
spike to the background profile p^, which is given either 
by the a-profile or by the Plummer profile, i.e. 



A_ 

7 



Ph(r-bh) 



(21) 



a condition which allows us to determine the constant 
A. Since the power law profile diverges at the centre, we 
must introduce a cutoff radius r m i n fixed either by the 
self-annihilation rate or by capture into the black hole. 
Using the former condition, the cutoff radius is given by 



A < a xx v > t hh lh 



(22) 



where ibh is the age of the SMBH, here taken to be equal 
to 10 10 yrs. For M31 the radius of the sphere of influence 
is about 7.7 pc. We have first considered a profile with a 
central core given by the Plummer profile for the back- 
ground, which results in a spike profile, p sp oc r -1 ' 5 . The 
7-ray flux was computed as before from (6), using (20) 
in the interval r m i n < r < rbh and the Plummer profile 
(14) for r > rbh- No significant enhancement is obtained 
whether the threshold energy is taken at 0.1 GeV or at 
1.0 GeV. However a different result is obtained if we use 
the a-profile (13) as the background profile. The spike 
profile is much steeper (7 ss 2.3), causing a more rapid 
increase of the density towards the centre and an impor- 
tant enhancement of the 7-ray flux. The strengthening 
of the signal is about a factor of 38 at 0.1 GeV and in- 
creases up to a factor of 55 at 1 GeV, due to the higher 
spatial resolution of GLAST at higher energies. As a con- 
sequence, if the halo of M31 has a profile with a central 
cusp, the 7-ray emission can be detected by GLAST at 
0.1 GeV if m x < 300 GeV and at 1.0 GeV if m x < 500 
GeV. These limits are considerably larger than those de- 
rived in Section V and open new possibilities to impose 
bounds on the neutralino mass and/or on the dark mat- 
ter distribution. 

The effects of a central spike are more dramatic in 
M87, since this galaxy has an extremely massive central 
black hole. The radius of influence of this SMBH is about 
97 pc. As we have seen in Section IV, the dark matter 
profile of M87 can be quite well represented by the a- 
profilc (13). Thus, the central spike has also a profile 
p sp oc r~ 2 3 , but the greater dark matter mass inside the 
sphere of influence produces an enhancement larger than 
that obtained for M31. The boost factor is about 200 at 
0.1 GeV and 234 at 1.0 GeV. These large factors lead to 
detectable signals either at 0.1 GeV if m x < 60 GeV or 
at 1.0 GeV if m x < 100 GeV. 



VIII. CONCLUSIONS 

Observational results favor a Universe whose main 
matter content is non-baryonic and dark. Presently, 
the most plausible DM candidate is neutralino; the 
lightest stable supersymmetric electrically neutral Majo- 
rana fcrmion. The possibility that neutralinos might be 
detected indirectly through their annihilation products 
such as 7-rays and neutrinos are explored by increasingly 
more sophisticated experiments which mainly search for 
energetic neutrinos from the centre of the Earth or the 
Sun and or for 7-rays from the galactic halo and various 
extragalactic sources. Although, the goal of DM detec- 
tion has not yet been achieved, upper bounds have been 
set on parameters such as neutralino mass. 

In this work, the 7-ray emission from neutralino anni- 
hilation has been studied. Our aim has been to use the 
results of indirect search experiments mainly GLAST to 
constrain the neutralino mass and the halo density pro- 
files as best as possible. Hence, we have used a simple de- 
scription of the annihilation process in which the only free 
parameter is the neutralino mass. We have estimated the 
photon production rate by fixing the present relic neu- 
tralino density to its equilibrium value at the decoupling 
temperature (f2 x = 0.26) favored by the latest CMB data 
and by assuming a QCD jet description of the annihila- 
tion process. Photon production rates above 0.1 and 1.0 
GeV are in the range 10~ 26 up to 10~ 24 ph cm 3 s _1 for the 
considered range of masses (10 < m x < 2000 GeV), con- 
sistent with results derived from supersymmetric codes 
such as DarkSUSY and SUSPECT. The mass interval 
that is explored here is compatible with the lower bound 
(<~ 30 GeV) given by the accelerator data and with the 
recent upper limit (~ 800 GeV) set by the elastic scat- 
tering of neutralinos from protons, including constraints 
imposed by new measurements of the anomalous mag- 
netic moment of the muon (g^-2) [69]. 

The DM density profile, was studied in the light of 
the different spatial resolution of the future experiment 
GLAST at different energies. Two different dark matter 
density profiles were adopted in our calculations: a profile 
with central cusp, in agreement with results of numerical 
simulations and one with a central core, more representa- 
tive of the observational data on rotation curves of bright 
spirals. These profiles are characterized by two parame- 
ters determined by the halo mass and by the mass density 
of dark matter at a given distance to the centre of the 
halo (in the case of the the Milky Way and M31) or us- 
ing data on the X-ray emissivity (in the case of M87). 
For dSphs, a modified Plummer profile, able to account 
for the observed raising velocity dispersion profile, was 
adopted for Draco whereas a King profile was assumed 
for Sagittarius. 

We have studied, by generating and experimenting 
with many different numerically simulated halos, the ef- 
fect of dumpiness on the 7-ray emission from M31 and 
conclude that the enhancement is less than a factor of 
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two for detectors centered on the object and having ac- 
ceptance solid angles of 10~ 5 — 10~ 3 sr _1 . The local en- 
hancement factor can be as large as 40-50 at the bound- 
aries of halos with central cusps. However, the absolute 
fluxes expected for a detector pointing away from the 
centre will be, in general, below the limit of detectabil- 
ity. These conclusions are in agreement with studies on 
the dumpiness effects either in the galactic halo [38,77] 
or in M31 [76], where maximum enhancement factors of 
about a factor of three were obtained. 

Our calculations show that if the dark matter halo 
of the Milky Way has a density profile with a central 
core, then the expect 7-ray intensity above 1 GeV in 
the direction b = 90° and averaged over a solid angle of 
0.5° x 0.5°, is comparable to the EGRET bounds [58], 
if m x < 50 GeV. However, a profile with a central cusp 
gives an expected intensity at least one order of magni- 
tude smaller, requiring an important enhancement due to 
the halo dumpiness to explain the data. This could be 
an additional argument in favor of high baryon-to-dark 
matter ratios in the central regions of bright galaxies, in 
agreement with analyses of rotation curves. 

Different spatial resolution at different energies of 
GLAST permits us to obtain informations either on the 
neutralino mass or on the dark matter distribution in 
the halos. Combining the informations on both energy 
thresholds, four possibilities arise which are written be- 
low. 

For M31, including effects of the spike produced by the 
central SMBH, the following situations may occur: 

a) If no detection is made at either energies then one 
may draw the conclusions that m x > 20 GeV and that 
the profile has a central core. 

b) If detection is made at both energy thresholds then 
the profile has a central cusp and m x < 300 GeV. 

c) If detection is made at 0.1 GeV but not at 1.0 GeV 
then we may conclude that the profile has a central core 
and m x < 20 GeV. 

d) Detection at 1.0 GeV but not at 0.1 GeV would result 
from a central cusp in the density profile and 300 < m x < 
500 GeV. 

A similar analysis for M87 yields: 

a) If no detection is made at cither energies then one may 
draw the conclusions that the "a-profile" is not valid at 
the centre where the density profile probably has a core 
or m x > 100 GeV. 

b) If detection is made at both energy thresholds then the 
"a-profile" is a good representation of the dark matter 
distribution near the centre and m x < 60 GeV. 

c) If detection is made at 0.1 GeV but not at 1.0 GeV 
then we may conclude that profile has a central core but 
no limits can be set on the neutralino mass unless an 
specific density profile is adopted. 

d) Detection at 1.0 GeV but not at 0.1 GeV would in- 
dicate that the "a-profile" is acceptable and 60 < m x < 
100 GeV. 

Because of their high dark matter content and their 
relative proximity the dSph Draco and Sagittarius are 



among the most favorable sources of 7-rays resulting from 
neutralino annihilation. However, we conclude that in 
spite of its large M/L ratio, the expected signals from 
Draco will not be detectable in the energy range of 
GLAST. On the contrary, our results suggest that an 
emission detectable by GLAST above 0.1 GeV is ex- 
pected for Sagittarius, if the neutralino mass is below 
50 GeV, a limit compatible with the LEP lower bound. 
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